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Abstract 

We describe a new type of top, the twisted top, obtained by appending a cocycle 
to the Lie-Poisson bracket for the charged heavy top, thus breaking its semidi- 
rect product structure. The twisted top has an integrable case that corresponds to 
the Lagrange (symmetric) top. We give a canonical description of the twisted top 
in terms of Euler angles. We also show by a numerical calculation of the largest 
Lyapunov exponent that the Kovalevskaya case of the twisted top is chaotic. 



1 Introduction 



We present a new top, called the twisted top, obtained by modifying the 
Lie-Poisson bracket for the charged heavy top. The charged heavy top, also 
introduced in this paper, is a heavy top [1-3] immersed in an electric field. The 
bracket for the charged heavy top arises from a semidirect product of SO (3) 
and M 3 xll The twisted top is not a top in the classical sense of a rigid body 
in a gravitational field. Rather, it is a mathematical construction obtained by 
using a different Lie group to build the Lie-Poisson bracket for the system. 
This abstract procedure is analogous to the manner in which tops are derived 
for SO(iV) [4], for SU(iV) (obtained in Hamiltonian truncations of the Euler 
equation [5,6]), and for other groups [7,8]. The construction method of the 
twisted top is also related to tops obtained by deformations of algebras [9]. 
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The bracket for the twisted top results from adding a cocycle to the charged 
heavy top bracket, so that the structure is no longer semidirect. Such brackets 
are classified in Thiffeault [10] and Thiffeault and Morrison [11], and the case 
we are considering is the simplest example of a Leibniz extension [12,10,11]. 
Because we are interested in how the nontrivial cocycle affects the dynamics 
of the system, we use the same Hamiltonian for the twisted top as for the 
heavy top. The bracket for the twisted top possesses three Casimir invariants, 
one of which differs from that possessed by the charged heavy top. 

A most interesting feature of the twisted top is that it retains the integrability 
property of the Lagrange top: it is integrable when it has an axis of symmetry 
(two moments of inertia are equal), its centre of rotation lies on the symmetry 
axis, and the electric field vanishes (or, equivalently, the top is uncharged). 
The conserved quantities are the energy, the angular momentum along the 
symmetry axis, and a third invariant which is a modification of the conserved 
component of the canonical momentum in the Lagrange case. 

The outline of this paper is as follows. In Section 2 we discuss the charged 
heavy top and describe its invariants and some of its integrable cases. In Sec- 
tion 3 we introduce the twisted top and its invariants. We show that it has an 
integrable case analogous to the Lagrange case of the heavy top. We give a 
canonical description of the twisted top in terms of Euler angles in Section 4. In 
canonical coordinates the difference between the twisted top and the charged 
heavy top is transferred from the bracket to the Hamiltonian, and appears as 
a term that can be interpreted as a momentum- dependent potential. Finally, 
in Section 5 we discuss our results and show by numerical calculation that the 
Kovalevskaya case of the twisted top is not integrable. 



2 The Charged Heavy Top 

Consider a heavy, charged top in constant gravitational and electric fields. 
The angular momentum vector is denoted by £, the position of the centre of 
mass is a vector a, and the position of the centre of charge is b. The direction 
and strength of the fixed gravitational and electric forces are given by the 
vectors a and (3, respectively. The frame of reference is the body frame, so 
that a and b are constant. The energy of such a top is 



where lo := I~ x t is the angular velocity and / is the moment of inertia tensor, 
which can be taken to be diagonal by an appropriate choice of frame. We 
assume that the top's rotation is slow enough that the magnetic fields set up 
by the motion of charges is negligible, and that the top is a perfect insulator, 



H(£, ct, (3) = \l ■ uj + a • a + (3 ■ b 
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so that the centre of charge remains fixed within the body. The charge, like 
the mass, does not have to be distributed uniformly, but only the centres of 
charge and mass couple to uniform gravitational and electric fields. 

The vectors a and f3, being fixed in space, rotate in the body frame. The 
dynamics of such a configuration can be generated by a Lie-Poisson bracket 
with a semidirect product structure, 

{/ > <?)sd =-^(VJxV £5 )-Q.(V £ /xV a9 + V a /xV £5 ) 

-/3-(VJxV^ + V/3/xV^), (2) 

where / and g are functions of (£,at,(3), and V is a gradient with respect to 
its subscript. Equation (2) is a simple extension of the bracket for the heavy 
top, which also has a semidirect product structure [2,3,13,14] (without a). 
The Casimir invariants of Eq. (2) are 

C 1 = H| 2 , C 2 = a-A C 3 =||/3|| 2 . 

The invariant C 2 says that the angle between a and j3 is constant, because 
by C\ and C3 their length is conserved. Therefore, the two vectors a. and (3 
fully describe the orientation of the rigid body, and there is a one-to-one 
mapping between ql and /3 and Euler angles. The phase space of the motion is 
thus SO(3) x IR 3 = T*SO(3). This is the same phase space as in the unreduced 
(canonical) system [15]. 

For the case with I x = I 2 , a = (0, 0, a 3 ) T , and b = 0, the charged heavy 
top reduces to the Lagrange top, also called the heavy symmetric top, and 
so is integrable (see for example Audin [16]). The invariants are the en- 
ergy H, £ 3 , and £ • ol. (There is also an integrable case with b = (0, 0, b^) T 
and (a x /3) • b = 0, i.e., where the forces and the centre of charge are copla- 
nar.) 



3 The Twisted Top 



In Thiffeault [10] and Thiffeault and Morrison [11], it is shown that the only 
bracket extension of two field variables (such as £ and a) is of the semidirect 
product type. To obtain an extension that is not semidirect, one requires at 
least three variables, which we take to be the same variables as for the charged 
heavy top. The simplest non-semidirect extension is then the Leibniz bracket 

{/ , <7} Le ib = {/ , 4)sd - e • ( V a / x V a g) , (3) 

where e is a parameter measuring the deviation from a semidirect bracket and 
is not necessarily small. Using the same Hamiltonian (1) as for the charged 
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heavy top in the bracket (3), we obtain the equations 



£ = {i,H\ 
a = {a. , H} 
(3 = {(3,H} 



Lcib 



Leib 



Lcib 



f xw + axa + /3xb, 
otxuj + eflxsi, 
(3 x u). 



(4) 
(5) 
(6) 



These are the equations for the twisted top. The term proportional to e in 
the 6t equation adds a "twist" which means that a does not simply rotate 
rigidly (though (3 still does). This is reflected in the Casimir invariants, which 
are now 



Since the length of a. is no longer preserved, the invariant C 2 does not im- 
ply that the angle between a. and f3 is constant. However, the length of the 
projection of a. onto j3 is preserved. 

For a positive-definite moment of inertia tensor, the energy surfaces of the 
twisted top are bounded, as can be seen from the following argument. First 
note that the components of £ cannot diverge without at or (3 also diverging, 
since £ enters the Hamiltonian in a positive-definite quadratic form. From the 
invariant C 3 , the components of (3 are finite. To have unbounded surfaces, and 
still conserve Ci, both ||£|| and ||a|| must go to infinity, with ||a|| 2 ~ —lei ■ f3. 
But with this functional relation it is not possible to have ||£|| — >• oo whilst 

1 1 2 

preserving the Hamiltonian H, since its kinetic part is proportional to \\£\\ and 

1 ' 1 1 1 /2 

its potential part to ||a|| ~ \\£\\ , precluding any balance. We conclude that 
the energy surfaces are bounded. This will be important in Section 5 where 
we try to demonstrate chaotic behaviour by computing the largest Lyapunov 
exponent. 

The twisted top also has an integrable Lagrange case. It is obtained, as for the 
charged heavy top, by letting I ± = I 2 , a = (0, 0, a 3 ) T , and b = 0. The energy H 
and the third component of the angular momentum l 3 are still conserved, 
whereas the third invariant becomes 



We call this integrable case the twisted Lagrange top. We can verify that P is 
conserved directly from the equations of motion (4)-(6), 

p = £ ■ a + £ ■ 6l + e ha 3 (3 3 

= (£ x uj) ■ ol + £ ■ (a x u) + e (3 x a) + e h a 3 ((3 x a;) 3 
= e a • {£ x (3) + e a • (J3 x h u>) = 0, 

where we equated JiWto£ in the last triple product because only the first two 
components of u> are involved, and Ii = I 2 - It is straightforward to verify that 



d= \\<x\\ 2 + 2e£-(3, C 2 = a-{3, C 3 = \\(3\\ 2 . 



(7) 



P = £ ■ a + e ha 3 (3 3 . 
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the invariants {H, £3, P} are in involution, i.e., they commute with respect to 
the bracket (3) — a necessary condition for integrability. The commutativity of 
the invariants carries over to the canonical variables of Section 4. 



4 Canonical Description 



Since the twisted top is a Hamiltonian system, there exists a coordinate trans- 
formation on the symplectice leaves (the constraint surfaces described by the 
Casimirs) that makes the system canonical. We now proceed to find such a 
coordinate transformation, in a manner analogous to the reduction of the rigid 
body and the heavy top [14,17]. The transformation we describe will be from 
the three Euler angles q= (4>,ifj,9) T and their corresponding canonical mo- 
menta p = (p<p,Pip,Pe) T (6 coordinates) to the vectors (£, a., (3) (9 coordinates, 
3 Casimirs). We show that the transformation is invertible on the symplectic 
leaves, so that it can be used to canonize the system. 

Following the heavy top reduction [14], since the vector (3 rotates rigidly 
(length conserved), it is fixed in the space frame, and we write 

(3 = A(q)w, (9) 

where the rotation matrix A is 

cos ip cos — cos 9 sin sin ip cos ip sin + cos 9 cos sin ip sin ip sin 9^ 
— sin ip cos — cos 9 sin cos ip — sin ip sin + cos 9 cos cos ip cos ip sin 9 ■ 

sin 9 cos cos 9 



sin 9 sin . 



The matrix A transforms vectors from the space frame to the body frame (we 
are following the convention of Goldstein [18, p. 147] for the definition of A.) 
The vector w is constant and fixed in space. Since rotations preserve lengths, 
we have C 3 = \\f3\\ 2 = w 2 , where w = ||w|| > 0. 



For the angular momentum, we take 

£ = L(q) p, 

where L is more concisely defined via its inverse, 



(10) 



L- 1 := 



^sin 9 sin tp sin 9 cos ip cos 9 ^ 
1 

cos ip — sin -0 



5 



This is the usual transformation one makes when reducing the rigid body [17, 
p. 499], where L maps T* SO (3) to 0*. 

Finally, for a we try the form 

a = A(q)v(q,p), (11) 

where we have allowed v to depend on the Euler angles and canonical momenta 
in an effort to conserve the Casimirs C\ and C 2 . We then have 

d = ct ■ (3 = v • A T A w = v ■ w, 

where we used the orthogonality of A. Decomposing v into a part v± e± per- 
pendicular to w and a part ve» parallel to w, we obtain d = vw. Since w is 
constant, we require v to also be constant. 

The norm of ct is 

||af = v • A T A v = v\ + v 2 = d - 2e£ ■ (3, 

where we used the definition (7) of d- We solve this for v\, 

v \ = d - v 2 - 2e£ ■ (3 = ||a|| 2 - ( "^f ^ > 0, (12) 

with v]_ — if and only if ex and (3 are collinear. For convenience, define the 
constant 



v := d - v = \\ct\\ - v ^' +2e£-(3, 

which we will use from now on instead of d- Then Eq. (11) becomes 

a = A vey + ^/r]-2e£- /3e ± (13) 

Note that the vectors o; and (3 are collinear (a x (3 — 0) if and only if rj = 2e I ■ (3. 
Assume that they are initially not collinear (77 7^ 2e £ ■ (3). The time evolution 
of 2e£ ■ (3 is obtained from (5) and the conservation of d, yielding 

2s j t {£- (3) = -j t ||a|| 2 = -2s a • (a x (3). (14) 

If rj = 2e £ ■ (3 initially, then it remains so for all times, because then the right- 
hand side of (14) vanishes. Conversely, if rj ^ 2e£- (3 initially, then the two 
vectors a and (3 are never collinear. 

This is crucial because it tells us that we can always invert Eqs. (9) and (13) 
for ((f), if;, 6), as long as a: and (3 are not initially collinear. The inversion is done 
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as follows: take w as the spatial z-axis. Then (3 is z', the transformed z-axis, 
which allows us to determine ip and 9, but not since it represents a rotation 
about the z-axis. We then use vj_ to define the x-axis, which allows us to find 
from a. (as long as 7^ 0, but we showed that it is sufficient to require this 
initially). But since the 2e£- (3 term only affects the magnitude of Vj_, not its 
orientation, we conclude that the Euler angles are only a function of oc and (3, 
not of £. We can then go back and solve (10) for the canonical momenta. 
(Provided det L~ l = sin 9 7^ 0, the coordinate singularity inherent to Euler 
angles. This singularity can be avoided by "inflating" the phase space [19].) 

It is straightforward to show that in the coordinates (q, p) the bracket (3) 
does indeed have the canonical form. 



In canonical coordinates, the 'potential' part of the Hamiltonian (1) becomes 
V(a,f3) = aa + (3h 



with 



The matrix L T A is 



(v a + w b) • A ey + J 7] - 2e £ ■ (3 a • A e± 



^•/3 = p-L T (q)A(q)w. 



^ - cot 6* sin cot 9 cos 1^ 



(15) 



L T A 



esc 6 sin — esc cos <j) 
cos sin (j) 0^ 



The integrable case of the twisted top has I\ = J 2 , a = (0, 0, a 3 ) T , b = 0, for 
which the kinetic energy is independent of (j) and -0, and the potential becomes 



V((f>, 6, p) = a 3 (v e 3 • A e,, + yj v - 2e £ ■ (3 e 3 • A e ± ) . 

Note that both L T A and e 3 • A = (sin 9 sin 0, — sin 6 cos </>, cos 6) are indepen- 
dent of -0, so that in the twisted Lagrange top case ip is cyclic (p^ conserved). 

A particularly simple choice is e\\ = e 3 , = e±, w = 1, for which the potential 
is 

V(4>,6,p) = a 3 (vcos9 + ^Jr) — 2e sin 9 sin cf>). 

Though simple, this description does not reduce nicely to the Lagrange top 
when s — > 0. The choice e_|_ = e 3 , ey — e±, w — 1, which gives 

V(<f>, 9, p) = a 3 (v sin 9 sin + \Jrj — 2e£ ■ (3 cos 9), 
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with 

t • (3 — —p^ cot 9 sin <fi + p^ esc 9 sin (/) + p e cos <p. 

has the Lagrange top form when v — e — 0, at the cost of a more complicated 
expression for £ • j3. 



5 Discussion 

We have introduced a simple generalisation of the heavy top by giving it 
charge and placing it in a constant electric field. By deforming the bracket of 
this charged heavy top, we have obtained a new top that we call twisted. We 
have found that the twisted top possesses an integrable case analogous to the 
Lagrange top. It is then natural to ask if such a deformation always preserves 
integrability. An affirmative answer would be very surprising, considering the 
delicate nature of integrable systems, and indeed it does not seem to be so for 
the system at hand. 

We investigate this by looking at the Kovalevskaya case of the twisted top. 
For the uncharged heavy top (b = 0, e = 0), the Kovalevskaya case involves 
setting Ii = I 2 , h = 2Ii, and a 3 = 0. This top is integrable [20,16]. The 
analogous case for the twisted top, as for the Lagrange top, simply involves 
a change of bracket by setting e ^ 0. Specifically, we choose e — 1 and a = 
(-1,0, Of. 

Figure 1 shows a plot of the instantaneous largest Lyapunov exponent of the 
twisted Kovalevskaya top, after averaging over 20, 000 random initial condi- 
tions integrated numerically (dotted line). The solid line is a least-squares fit to 
help determine the Lyapunov exponent to greater accuracy, using an asymp- 
totic form of the averaged exponent [21]. The Lyapunov exponent is A ~ 0.122, 
suggesting that the twisted Kovalevskaya top is chaotic. For comparison, 
the dashed line shows the same calculation for the untwisted (ordinary) Ko- 
valevskaya top, which shows the Lyapunov exponent going to zero. Thus, the 
twisted case appears to be chaotic, whilst the untwisted case is not. We con- 
clude that integrability does not always survive deformation, contrary to the 
Lagrange case. Note that we can infer chaos from a positive Lyapunov expo- 
nent because we showed in Section 3 that the motion takes place in a bounded 
region of phase space [22]. The presence of chaos does not rule out the exis- 
tence of an integrable case with parameters "close" to the Kovalevskaya values 
(in the sense of differing only by terms involving e). We have not found such 
a case. Another approach would be to use a Kovalevskaya-Painleve analysis 
to determine other integrable limits of the twisted top. 

A rigorous demonstration that the twisted Kovalevskaya top is chaotic could 
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in principle be achieved using a Melnikov analysis, as was done by Holmes 
and Marsden [14] for the heavy top. The present case is less straightforward 
because of the complicated form of the homoclinic orbits. 

It would of course be of great value to find a physical realisation of the 
twisted top. One could use either the noncanonical picture, Eqs (4)-(6), or 
the canonical picture, given by the standard free rigid body Hamiltonian [14] 
with Eq. (15) for a potential. Regardless of the physical interpretation, the 
twisted top remains an object worthy of study in its own right, because of its 
interesting integrable case and peculiar geometry. It would be worthwhile to 
carry out a topological classification of the bifurcations of the phases space of 
twisted Lagrange top, as was done by Dullin et al. [23] for the Kovalevskaya 
top. This is complicated by the need to find a good surface to make Poincare 
sections in the canonical coordinate space. 
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Fig. 1. Lyapunov exponent for the twisted top, averaged over initial conditions. 
The dotted line is for the Kovalevskaya top and the solid line is the func- 
tion 0.497/t + 0.113/v^ + 0.122, obtained by a least -squares fit and yielding the 
value A ~ 0.122 (See Ref. [21]). For comparison, the dashed line is the averaged 
Lyapunov exponent for the ordinary Kovalevskaya top. 
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